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ABSTRACT 

Fisher matrices are used frequently in the analysis of combining cosmological constraints from 
various data sets. They encode the Gaussian uncertainties of multiple variables. They are simple to 
use, and I show how to get up and running with them quickly. Python software is also provided. I 
cover how to obtain confidence ellipses, add data sets, apply priors, marginalize, transform variables, 
and even calculate your own Fisher matrices. This treatment is not new, but I aim to provide a clear 
and concise reference guide. I also provide references and links to more sophisticated treatments and 
software. 

Subject headings: cosmology 



1. OUTLINE 

I explain how to do/obtain the following with/from 
Fisher matrices: 

§ [2] Confidence Ellipses 
§ jsj Manipulation: 

Marginalization, Priors, Adding Data Sets 

4 How to Calculate your Own Fisher Matrices 

5 How to transform variables 

6 Dark energy pivot redshift 

7 Discussion (brief) about what Fisher matrices are 
5 Software I've come across (including my own) 
^ How you can contribute to this paper 

2. FISHER MATRICES => CONFIDENCE ELLIPSES 

The inverse of the Fisher matrix is the covariance 
matrix: 



\F]-' = [C] 



xy 



(1) 



(Ta; and ay are the 1-a uncertainties in your parameters x 
and y, respectively (marginalizing over the other), a. 



xy 



pax^y, where p is known as the correlation coefficient, p 
varies from (independent) to 1 (completely correlated). 
Examples are plotted in Fig. [l] 
The ellipse p arameters are calculated as follows (e.g.. 



Unknown 



2008): 
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We then multiply the axis lengths a and & by a coefficient 
a depending on the confidence level we are interested in. 
For 68.3% CL (I-ct), Ax^ « 2.3, a = y/A^ « 1.52. 
Other values can be fo und in Table [l] These can be 
calculated following e.g., Lampton et al. ( |1976 1. 
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TABLE 1 
Confidence Ellipses: 



a CL 


Ax^ 


a 


1 68.3% 

2 95.4% 

3 99.7% 


2.3 
6.17 
11.8 


1.52 
2.48 

3.44 


The area of the ellipse is given 


by 




A = TT(aa){ab) 
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= 7r(Ax^)a5 




(6) 




-p2 


(7) 



The inverse of the area is a good measure of fig ure of 
merit The Dark Energy Task Force (DETF; Alb rechtj 
et al.||2006| |2009|) used FOM = tt/A for the abilit"^^ 



experiments (WL, SN, BAO, CL) to constrain the dark 
energy equation of state parameters {wq, Wa)- 

2.1. Probability P{x,y) 

Interested in the probability that specific values 
are correct for parameters x and yl The probability 
function P{x,y) given best fit values (a;o,?/o) and 1-a 
uncertainties {ax,ay) is calculated as follows: 
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l-p2 
P{x, y) = exp 
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with Aa; = x — Xq and Ay 
(uncorrelated x and y), the x 
For correlated x and y (p > 0), x 



y ~ yo- Note for p = 
2 formula looks familiar. 
^ is reduced. 



3. MANIPULATION: MARGINALIZATION, PRIORS, 
ADDING DATA SETS, AND MORE 

Consider a Fisher matrix provided by the DETF (Table 
[2]) for optimistic Stage IV BAO observations for the fol- 
lowing variables: {uJm,^A,^k), where ujrn = ^mh^ and 



Fig. 1. — 68.3% (l-c) confidence ellipses for parameters x and y with 1-cr uncertainties and ay and correlation coefl&cient p. In the 
first three panels, we plot as dashed lines the marginalized 1-a uncertainty for each variable: aa^ and aay, where a a \/2.3 ~ 1.52. In the 
bottom-right panel, we zoom in to show the intersections with the axes: zt/Jo-j, and it/3cry, where /3 2.13v'l — p (for p 1). 



TABLE 2 
Example Fisher Matrix 













2,376,145 


796,031 


615,114 




796,031 


274,627 


217,371 




615,114 


217,371 


178,014 






TABLE 3 






Corresponding Covariance Matrix 














3.20E-5 


-1.56E-4 


8.02E-5 


Ha 


-1.56E-4 


8.71E-4 


-5.25E-4 




8.02E-5 


-5.25E-4 


3.69E-4 



fl„i + ^A + ^k = 1- The covariance matrix (inverse of the 
Fisher matrix) is given in Table |3] For example, the top- 
left element tells us that Atj,„ « 0.00566 « V3.20E~5. 



3.1. Marginalization 

When quoting these uncertainties on lu^, the other 
variables (f^Ajf^fc) have automatically been marginalized 
over. That is, their probabilities have been integrated 
over: they have been set free to hold any values while we 
calculate the range of acceptable uj,n- 

To calculate a new Fisher matrix marginalized over any 



TABLE 4 
Fisher Matrix with Fixed Qj. = 







Ua 




2,376,145 


796,031 


Qa 


796,031 


274,627 



variable, simply remove that variable's row and column 
from the covariance matrix, and take the inverse of that 
to yield the new Fisher matrix. 



3.2. Fixing Parameters 

Suppose instead want the opposite: perfect knowledge 
of a parameter. For example, we want to consider a 
flat universe with a fixed value of fl^ = 0. To do this, 
simply remove ilk from the Fisher matrix (Table |4j). The 
new covariance matrix and parameter uncertainties are 
calculated from the revised Fisher matrix. 

Alternatively, the on-diagonal element corresponding 
to that parameter can be set to a very large value. For 
example, if we set the bottom-right element in Table |2] 
to 10^^, that would correspond to a 10^^ uncertainty in 
LOm, or nearly fixed. Note that higher values in the Fisher 
matrix correspond to higher certainty. 

3.3. Priors 

Rather than fixing a parameter to an exact value, we 
may want to place a prior such as Aflk = 0.01 (l-c). 
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In this case, simply add 1/<t^ = 10* to the on-diagonal 
element corresponding to that variable (in this case, the 
bottom left element). 

3.4. Adding Data Sets 

To combine constraints from multiple experiments, 
simply add their Fisher matrices: F = Fi + F2. Strictly 
speaking, any marginalization should be performed af- 
ter the addition. But if the "nuisance parameters" are 
uncorrelated between the two data sets, then marginal- 
ization may be performed before the addition. 

4. HOW TO CALCULATE YOUR OWN FISHER MATRICES 

Given the badness of fit x^i^iU): your 2-D Fisher 
matrix can be calculated as follows: 



[F] = 



1 



dxdy 



dxdy 



(10) 



In other words, = 2 
These derivatives are simple to calculate numerically: 

Q2^2 ^2(^^ ^_ y^-^ _ 2x^{xq, yo) + x^{xa - Aa;, yo) 



dx'^ 



dx 



X^ixo + Ax, yo) - x'^ixo - Ax, yo) 



2Aa 



d^X^ 



d 



dx^ 
dx 



dxdy dy 
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5. TRANSFORMATION OF VARIABLES 

Suppose we are given a Fisher matrix in terms of 
variables p — {x,y, z) but we are interested in constraints 
on related variables p' — (a, 6, c). We can obtain a new 
Fisher matrix as follows: 



where Mij = 

dp'j 
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and [Af]^ is the transpose. 

All of these partial derivatives should be evaluated nu- 
merically, plugging in best-fit values of the parameters. 

5.1. Transformation Example 

Suppose we are given a Fisher matrix in terms of 
(wm, r^Ai ^k)i but we are interested in (f2m, h). Here 
LUm = flmh^ and flk = 1 — — Suppose further that 
the best-fit cosmology is {Qjn,^A,h.) = (0.3,0.7,0.7). 
Our transformation matrix is evaluated as follows: 

duJm doJm duJri 



dilm dflA dh 
di^A dflA d^A 
dilm di^A dh 
dQ.k dQ.k dQ.k 
d^lm dilA dh 



2n^h 
1 

-1-1 



0.49 0.42 
1 

-1 -1 



(20) 

(21) 
(22) 



6. PIVOT REDSHIFT 

Given the dark energy equation of state parameteriza- 
tion 



w = Wo + (1 — a)wa 



(23) 



F' 



dpi dpj 



^ dp',n dp'n 



F,, 
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Let's spell this out explicitly. Here is the expression for 
element (a, b) in the new Fisher matrix: 



''"'da 96 5a 96 



dydXp 
da db 



9y p 
da db 



dz dx 



dz dy 



F,. 



dx dz 

da db 
dy dz 

dadb"^' 
dz dz ^ 

zz 



da db ^"^ ' da db da db 

This can be calculated using matrices: 
[F']^[Mr[F][M] 



^ XZ 

Fy 
F, 



(15) 
(16) 
(17) 

(18) 



where 1/a = 1 + z, if you have calculated a Fisher 
Matrix for dark energy parameters wq and Wa , go ahead 
and calculate the pivot redshift, too: 



-1 



1 



AWg 

pAwa 



(24) 



A t this redshift, w{z) is b est constrained (e.g.. Fig. 16 
of Huterer & Turner||2001 1. Rather than presenting con- 



straints on {'Wo,'Wa), constraints on {wp,Wa) can be pre- 
sented. That is, we constrain the value oi w aX z = Zp 
rather than at z = (along with w's rate of change with 
time Wa)- 

The {wp,Wa) confidence ellipse has no tilt; there 
is no correlation between the two, by definition.^ 

^ Thus the DETF chooses a more interesting elUpse to plot: 
{wp, Qde)- 
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But the area of the (wpjWa) eUipse is equal to the area 
of the {wo, Wa) eUipse. From this and Eq.^it follows that 



(25) 



And if w is constant, then Awp = I^wq. 



Derivation of the p ivot redshift formula follows from 
( Albrecht et al.|[2006[ ), calculating the uncertainty of 



Vp = Wq + {I - ap)Wa 



{Awp)^ = {Awof + ((1 - ap)Awaf + 2(1 



(26) 



ap)Awo,a 
(27) 



where Awo.a — pAwoAwa, and then minimizing Awp for 



7. DISCUSSION 

Fisher matrices encode the Gaussian uncertainties in a 
number of parameters. Confidence ellipses can be easily 
calculated over any pair of parameters. These provide an 
optimistic approximation to the true probability distri- 
bution. The true uncertainties may be larger and non- 
Gaussian. Note the best fit values themselves are not 
encoded in the Fisher matrices, and must be provided 
separately. 

Fisher matrices allow one to easily manipulate param- 
eter constraints over many variables. It is easy to add 
data sets, add priors, marginalize over parameters, and 
transform variables, as shown here. 

A more in-depth discussion of Fisher matric es and is- 
sues surrou nding their use can be found in (Albrecht 



et al. 



2009 



'i'his is the paper I'd wished I could find when I began 
my work with Fisher matrices: projections for cosm olog- 
ical constraints from gravitational lens time delays ( Coe 
(fc Moustakas||2609| ). 



http : // www, its . caltech . edu/'/.TEcoe/Fi 
http : / A?ww. physics .ucdavis . ed u/DETF. 



http : //www. cosmology . org. za/ 



isherk 
ast/f 



8. SOFTWARE 

Fisher. py^ Python - simple manipulation of Fisher 
matrices and p lotting of ellipse s 

DETFast^ ^ 
expectations 



|Albrecht et al.||2006[ ) JAVA - Compare 
of cosmological constraints from different 
experiments with your choice of pr iors with a few clicks! 



Fisher4Cast'' ( |Bassett et al.|2009^ Matlab - most so- 
phisticated 
Your ad here. 



9. CONTRIBUTE 

This is meant to be a brief guide, but if I've failed to 
reference another useful guide or your software or if I've 
neglected some detail (subtle or otherwise) about Fisher 
matrices, please e-mail me at coe(at)caltech.edu, and I'll 
be happy to update this document. Also please tell me 
if any section is unclear. 

If I have not covered a useful topic, it is probably 
outside my knowledge of Fisher matrices. For exam- 
ple, I have not covered the analysis of Monte Carlo 
Markov Chains (MCMC) as provided, for example, by 
the WMAP Lambda website.'^ If a generous reader could 
explain to me (or point me to an appropriate reference 
on) how to extract confidence contours and a Fisher ma- 
trix from a MCMC, I would be grateful and include the 
explanation here, giving due credit to the contributor. 



I thank Olivier Dore for referring me to the DETFast 
software written by Jason Dick and Lloyd Knox whom I 
also thank for answering my questions about their soft- 
ware. It is a valuable resource. Once I took off these 
training wheels and began to produce my own plots, 
DETFast is still a valuable resource for Fisher matrices 
calculated by the DETF encoding their estimates of cos- 
mological constraints from various future experiments. 

This work was carried out at Jet Propulsion Labora- 
tory, California Institute of Technology, under a contract 
with NASA. 



^ http://lambda.gsfc.nasa.gov/ 
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